Orientational pinning of quantum Hall striped phase 

E. E. TakhtamirovS and V. A. Volkovi 
Institute of Radioengineering and Electronics of RAS, 
11 Mokhovaya, Moscow 103907, Russia 



In ultra-clean 2D electron systems on (001) GaAs/AlGaAs 
upon filling high Landau levels, it was recently observed a new 
class of collective states, which can be related to the sponta- 
neous formation of a charge density wave ("striped phase"). 
We address to the following unsolved problem: what is the 
reason for stripe pinning along the crystallographic direction 
[110]? It is shown that in a single heterojunction (001) A3B5 
the effective mass of 2D electrons is anisotropic. This natu- 
ral anisotropy is due to the reduced (C2v) symmetry of the 
heterojunction and, even being weak (0.1%), can govern the 
stripe direction. A magnetic field parallel to the interface in- 
duces "magnetic" anisotropy of the effective mass. The com- 
petition of these two types of anisotropy provides quantitative 
description of the experiment. 



I. INTRODUCTION 

It had been assumed (l) in 1979 that a uniform 2D 
electron system in strong magnetic fields corresponding 
to the filling of the lowest Landau level, N = (filling 
factor v < 1), can be unstable against the formation of 
a ID charge density wave with a period of the order of 
magnetic length. This instability is due to the exchange 
interaction resulting in the effective attraction between 
electrons. The analysis was carried in the Hartree-Fock 
approximation, which overestimates the exchange inter- 
action and ignores electron-electron correlations. After 
the discovery of the fractional quantum Hall effect, it 
became clear that, for v < 1, it is the correlation inter- 
action that leads to the formation of a uniform state of 
the Laughlin liquid type. Nevertheless, the role of corre- 
lations diminishes with filling a large number of Landau 
levels and, in principle, one can expect the indicated in- 
stability to appear. It was predicted in 1996 ^ that a 
ID charge density wave may appear near the half-filling 
of Landau levels, beginning with N = 2 [pi. Such a 
striped phase with a period of order of the Larmor di- 
ameter should be energetically more favorable than the 
Laughlin liquid and the Wigner crystal j|||,0] . 

How should this phase be manifested in the transport 
measurements if it is really formed and pinned for some 
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reason? Similar problem was considered, probably, for 
the first time in a series of old papers j5UqM| , where the 
anisotropic conductivity was calculated |5|] and the high- 
frequency H and heating Q] effects were studied. In the 
presence of a periodic ID potential U(x) induced by the 
charge density wave, each Landau level transforms into 
a narrow ID band. At the edges of this band, the den- 
sity of states has a power divergence, which is cut off 
upon including weak scattering. As a result, the den- 
sity of states S(Ef) at the Fermi level has the shape of a 
two-teeth fork: a minimum in the center of the band and 
two peaks at the band edges, the lower the density rii m of 
scatterers, the higher the peaks. Incomplete filling of this 
band results in the formation of stripes differing in the v 
value [of the type vjiv — l)/v/ ..] and aligned with the 
y-axis. The transverse (a xx ) and longitudinal (o yy ) con- 
ductivities obey different mechanisms and qualitatively 
differently depend on v: the a yy conductivity is high, 
has the band character, and is inversely proportional to 
rii m and S 2 (Ef), while a xx is low, has the hopping char- 
acter, and is proportional to n im and S 2 (Ef), with the 
product a xx a yy being independent of scattering. These 
results were in fact confirmed and generalized in [pL[j]|. 
What did the experiment actually reveal? 

In 1999, the conductivity of an electron system with 
ultrahigh mobility in the (001) GaAs/AlGaAs structures 
was studied at very low temperatures near half-integer 
v > 9/2, and a new state was revealed ]To| , [TT| , |T^ , p~3| | 
which was assumed to be just the one associated with 
the formation of the striped phase predicted in This 
assumption was primarily based on the observation of a 
giant resistance anisotropy in this system. The ratio of 
resistances along the crystallographic directions [110] and 
[110] reaches the values of R xx /R yy ~ 5-3500, depend- 
ing on the sample geometry, where [110] is the "easy" 
conductivity direction. Moreover, the behavior of all 
conductivity tensor components qualitatively agrees with 
the theory j|||f|. The predicted behavior of the a xx a yy 
product near the half filling of the upper Landau level 
numerically agrees with the experiment It was also 
shown in Jl2|,[l3| that the magnetic field £?|| ~ 1 T parallel 
to the interface can change the direction of easy conduc- 
tivity. The authors of (l^] concluded that, at high enough 
J5|| , the direction of easy conductivity is perpendicular to 
the -By direction. Similar result was obtained in |l3f for 
B || || [110] near all half-integer v > 9/2 and for B"llO] 
near v — 11/2 and v = 15/2. The theoretical analysis 
|l5],[L6| of the influence of B|| carried out in the Hartree- 
Fock approximation in the model of parabolic quantum 
well partially explained the results. All this is strong ev- 
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idence for the formation of a striped phase. However, a 
mechanism responsible for the orientation of the charge 
density stripes along a certain preferred direction in a 
macroscopic sample (orientational pinning) at £?| = 
remains to be clarified. This is one of the fundamental 
unsolved problems. 

In this work, the Kroemer's assumption jL7) that the 
reduced (C^) symmetry of the potential in a heterostruc- 
ture based on semiconductors without inversion center 
can cause an appearance of a preferred direction for the 
conductivity is justified. Symmetry reduction means that 
the cubic axis normal to the interface is transformed from 
the fourfold mirror-rotational axis (54) to the twofold 
axis (C2). We demonstrate below that, owing to the 
asymmetry of the potential of an atomically sharp het- 
erojunction, the effective mass (EM) of 2D electrons is 
anisotropic (natural anisotropy). At the same time, the 
presence of By also gives rise to the EM anisotropy (mag- 
netic anisotropy) |18| . Therefore, the results of many- 
particle numerical calculations of the B|| effect on the 
orientation of striped phase Jjl],[l(| is natural to treat (to 
the lowest order in B?A as a manifestation of the mag- 
netic anisotropy of EM. Thus, the many-particle problem 
of orientational pinning of the striped phase reduces to a 
one-particle problem of determining the EM anisotropy. 
We derive below analytical expressions for both types of 
EM anisotropy (natural and magnetic) and demonstrate 
that they can compete with each other. At a certain mag- 
nitude and direction of By , these two types of anisotropy 
exactly cancel, leading to the disappearance of resistance 
anisotropy, in agreement with the experiment. 



II. NATURAL ANISOTROPY 

Before proceeding to the many-particle problem, it is 
necessary to obtain one-particle Hamiltonian for the con- 
duction band in a (001) A3B5 heterostructure. As was 
shown in |H| , the correctly constructed multiband set of 
equations for the envelope functions retains information 
on the heterostructure symmetry (C^), which is lower 
than the symmetry Td of the constituent materials. This 
symmetry reduction is described, in particular, by cer- 
tain short-range potentials localized at the heterointcr- 
face. Mixing of heavy and light holes at the center of 
the 2D Brillouin zone is one of the consequences of sym- 
metry reduction []l9| , po|| . This mixing explains giant op- 
tical anisotropy (with the same principal axes [110] and 
[110]) that was discovered in |2l]] for the quantum wells 
based on semiconductors with different cations and an- 
ions. Evidently, low symmetry should also manifest itself 
in the equation for the envelope functions in the con- 
duction band. Nevertheless, a single-band equation ob- 
tained in [ 19j carries no information on the C% v symmetry 
because the corresponding small contributions were ne- 
glected. We must now take them into account. Since 
the terms with the symmetry higher than Ci v are of no 



interest here, the effective Hamiltonian can include only 
the operators of kinetic and potential energies used in the 
standard EM approximation, as well as the anisotropic 
contribution of C% v symmetry, which will be obtained 
below. 

A single-band Hamiltonian of the Cn v symmetry can 
be obtained by the method of invariants. Leading aside 
spin-orbit interaction, one can conclude that the C% v 
symmetry should manifest itself in the kinetic energy op- 
erator. Let us direct the 2D quasi-momentum compo- 
nents along the cubic axes: j?^ || [100] and p y \\ [010] (z axis 
is along [001]). Then the part of kinetic energy opera- 
tor quadratic in the generalized 2D momentum (P x , P y ) 
should be 
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Here, m* is the EM of conduction band, and the quantity 
A (which may depend on z) accounts for the natural 
anisotropy of EM in the plane of 2D electron gas. Let 
us obtain explicit expression for A by using a multiband 
matrix Hamiltonian [ ft9| acting on the column of envelope 
functions. In the k representation, it takes the form 
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where n and n! are the band indices. The first term in 
Eq. (||) includes the contributions from the smooth po- 
tentials and kp interaction and has the standard form. 
The second term in Eq. (|^) is a contribution from the 
atomically sharp heterointerface potential taken to a first 
order in the small parameter ka, where l/k is the char- 
acteristic length of changing the envelope functions and 
a is the lattice constant. One can pass to the single- 
band variant of envelope-function method by applying 
the perturbation theory, with the kp and Do operators 
as perturbation. The second order in the kp interaction 
gives the first (standard) term in Eq. (|l|). The third order 
(second order in kp and first in D ) provides the second 
term of Eq. (|l|), with A(z) = aS(z) and 
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Here 5{z) is the Dirac 5-function, z = defines the het- 
erointerface position, (n \ pi \ n') is the ith component of 
the interband momentum matrix element, c is the index 
of conduction band, mo is the mass of free electron, and 
e„ is the energy of the nth band edge in one of the struc- 
ture materials. In the simplest model, the key parameters 
of the theory (£>onn') have the form 
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The functions G(z) and SU(r) are so defined that the 
crystal potential of the heterostructure has the form 
U(r) = Ui{r) + G(z)SU(r), where Ui and U 2 = Ui + 5U 
are the crystal potentials of the structure materials. Note 
that the Doxy parameter accounts for the mixing of the 
heavy and light holes at the center of the 2D Brillouin 
zone, with X and Y being the indices of the Bloch func- 
tions corresponding to the edge of the T 15 valence band 
and transforming as x and y under symmetry operations 
of the T d group MM. 



III. INCLUSION OF THE MAGNETIC 
ANISOTROPY 



Reducing the tensor of reciprocal effective mass to the 
principal axes, so that a; || [110] and j/ 1 1 [110] in the new co- 
ordinates, and introducing magnetic field B in the vector- 
potential gauge A = (B y z, —B x z + B z x,0), one obtains 
for the orbital part of the 3D Hamiltonian of conduction 
band 
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Here, V(z) is the effective potential of the conduction 
band edge, e is the elementary charge, and c is the light 
speed (we hope, it will be no confusion with the conduc- 
tion band index). For the finite thickness of the 2D layer, 
the magnetic field component parallel to the heterointer- 
face can be treated perturbatively [0. To second order 
in B\\ , this procedure brings about a diamagnetic shift of 
the dimensional-quantization subbands and an increase 
(for the lowest subband) in EM in the direction perpen- 
dicular to By. The natural EM anisotropy also can be 
treated perturbatively. For simplicity, we assume that 
B|| is parallel to either [110] or [110], so that B x B y = 0. 
Collecting all terms second-order in B^ and first-order in 
a, one obtains the following expression for the 2D Hamil- 
tonian of the lowest subband: 



The parameters of the natural EM anisotropy and the 
EM anisotropy induced by the magnetic field are 
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where E m is the energy of the bottom of the mth subband 
at B — 0. The expression for in Eq. ^ is valid 
to the terms second-order in the parameter hu c /{E2 — 
Ei), where uj c = eB z /m*c. For the field B z — 2.5 T (in 
the experiment jl3| this field corresponds to the filling 
factor v = 9/2), one can neglect this correction in the 
estimation of A^, because hoj c w 4meV, while the gap 
E2 — Ei should exceed the Fermi energy Ef measured 
from the lowest subband; one has Ef w lOmeV for the 
2D electron concentration N„ — 2.7 x 10 11 cm~ 2 . 



IV. ESTIMATES 

Based on the experimental data |l3|] , we estimate A Ilat 
and A B for B\\ = 0.5 T (if B||||[110], this magnetic field 
converts the resistance from anisotropic to isotropic; at 
larger B» the direction of "easy" conductivity rotates by 
90°). Since the information on the samples is incom- 
plete, we carried out a series of self-consistent calcula- 
tions by varying the concentration N a of residual accep- 
tors in GaAs. At iV a 
]l3[ , the results are 
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For other N a values (from 10 13 to 10 15 cm -3 ), the results 
differ from Eq. (||) and (||) by a factor less than two. 
Equations (^) and (||) yield the following value for the 
EM anisotropy induced by magnetic field Bu — 0.5 T: 



A B = 1.3 x 10 -3 = 0.13%. 



(10) 



The parameter a can be determined from the equality 
A nat = A B to give 

0.65 xlO- 8 cm 1n20 
a = =1.1x10 cm/j. (11) 
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The two-band approximation with energy gap E g yields 
the following estimate for Eq. (||): 
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Thus, it follows from the experimental data ||13(] and 
Eqs. @ and @ that Doxy ~ 0.4 x 10" 8 eVcm. Let 
us compare this value with the literature data. 

It was found in [^0| that different estimates carried out 
for the GaAs/AlAs heterostructures either on the basis 
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of pseudopotential and tight-binding calculations or from 
the comparison with the experiment on the anisotropic 
exchange splitting of exciton levels in the II-type super- 
lattices GaAs/AlAs lead to a sizable scatter of the Doxy 
parameter. The value obtained in p0| lies in the range 
(0.35, 0.99) x 10 -8 eVcm. Using linear interpolation, one 
obtains the upper bound Doxy = 0.3 x 10 _8 eVcm for 
the GaAsZAlo.3Gao.7As heterostructure. This value is in 
a fair agreement with the value obtained above. 



V. DISCUSSION 

We can now conclude that the natural anisotropy of 
EM is likely the mechanism that pins the stripe directions 
at B\\ = (see also the end of section 1). It follows from 
this conclusion that the parameter a entering Eq. (Q) is 
negative, a < 0. The competition between the natural 
anisotropy A nat and anisotropy induced by the mag- 
netic field By = (0, By) makes the 2D electron spectrum 
at Bn = 0.5 T isotropic. As a result, the stripe directions 
are randomized and the resistance becomes isotropic. On 
further increase in By , the magnetic anisotropy prevails 
and the stripes rotate at 90° . The role of EM anisotropy 
in the formation of many-electron anisotropic states can 
be understood as follows. A 2D electron system with 
anisotropic EM and isotropic Coulomb interaction is, ob- 
viously, equivalent to a 2D electron system with isotropic 
(cyclotron) mass and anisotropic Coulomb interaction. 
One can expect that this effective anisotropic interac- 
tion is precisely the one which pins the orientation of the 
striped phase to ensure its observation in the magneto- 
transport. 

For holes, the heterointerface contribution of symme- 
try C2v (and, hence, responsible for the pinning of the 
striped phase) is greater than for the conduction band, 
because it appears in the first-order perturbation treat- 
ment ]l9| ), whereas the anisotropic EM in Eq. (|l]) was 
obtained in the third order. For this reason, one would 
expect that the hole striped phase is more stable and can 
form upon filling the lower Landau levels (cf. [0). 
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